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Project Concept 


■ AIRS radiances currently assimilated operationally in GFS and NAM 

• Cloud-free radiances from 281 -channel subset 

• Cloud checks performed within GSI to determine which channels peak above cloud top 

• Inaccuracies may lead to less radiances assimilated or introduction of biases in cloud- 
contaminated radiances 


■ Use AIRS L2 retrieved profiles to better understand the optimal three-dimensional 
distribution of AIRS radiances assimilated within GSI to engage the operational 
DA community regarding strategies for assimilating hyperspectral radiances 


• Cloud contamination, channel reduction, spatial data reduction 
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■ Lowest extent of quality 
AIRS L2 profiles determined 
by quality indicators in 
preprocessing 

■ Use MODIS as an additional 
resource to determine cloud 
location and vertical extent 


Lowest extent of useable L2 profiles 
determined in preprocessing 


Lowest extent of useable LI B 
radiances determined within GSI 


5 


Experimental Setup 

Developmental Testbed Center (DTC) GSIv3.0 and WRF-NMMv3.3 code 
configured in forecast cycling methodology that mimics the operational NAM 

Real-time BUFR files archived during assimilation period (4 Nov. -20 Dec. 2011) 

• Satellite: AIRS, AMSU, HIRS, MHS, GOES Sounder, GPSRO, radar winds 

• Conventional: All observations used in EMC's Table 4 

Two "parallel" 4-week experiments with 2-week spin-up: 


RAD 

o assimilate AIRS radiance data 
using operational procedures 

PROF: 

o append PREPBUFR to include 
AIRS profiles as sondes 
ensuring consistency with 
real-time RAD swath locations 

o quality flag P best to select data 
in the vertical to be assimilated 

o no observation thinning 


Schematic for GSI scripts (DiMego, personal communication, 2011) 
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Bulk Cloud Information 


Mean MODIS cloud fraction: 20 Nov - 20 Dec 2011 




Mean MODIS CTP: 20 Nov - 20 Dec 2011 


High cloud tops 


low cloud tops 


■ Mean cloud information from Aqua MODIS interpolated to WRF-NMM grid 

■ Main focus is on strip of persistent low, opaque clouds just north of Equator 

■ Regions with low, opaque clouds (assimilate additional channels above cloud) 
and regions with cloud gradients (assimilate in partly cloudy or scene incorrectly 
deemed cloudy) should be areas where profiles will be most impactful 
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Forecast Impact 


Mean (21 Nov-19 Dec 2011) 500 hPa ZAC difference Mean (21 Nov-19 Dec 2011) 500 hPa T AC difference 
at F48 for 00Z initialization at F48 for 00Z initialization 



Region of >0.10 difference 
along equator downstream of 
low cloud region 


■ Using same-cycle analysis valid at forecast time as verification field 


■ NCEP GFS climatology interpolated to NMM grid used for AC calculation 


500 hPa anomaly correlation differences between profile and radiance — 

Downstream of low cloud region, which is coincident observations valid for the 00 
and 12 UTC analyses 
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Impact Difference for Select Case 



- 1.8 - 1.5 - 1 2 - 0.9 


- 0.3 0.3 0.6 0.9 1.2 1.5 1.8 


Temperature (K) ID at cr=39 (~500 hPa) for 00Z analysis on 22 

November 2011 


■ Impact Difference (ID) was 
calculated for each 00Z analysis 
and interesting cases for 
further investigation were 
selected 

ID t j = | RADALYSi j - RADBKGDi j\ 

- \PROFALYSi j - PROFBKGDi 

■ What follows is an example of 
the analysis being performed 
for a single case (22 Nov 2011) 

■ Following slides examine 
possible explanations in GSI 
diagnostics and MODIS cloud 
products for area over SE 
Pacific near the equator to help 
explain improved profile 
forecasts 
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Comparison to MODIS CTP 

■ Two regions with «1.5K larger analysis impact in profile analysis 

■ Overall, GSI does a good job of determining cloud top pressure (CTP); devil is in 
the details 

■ For regions of largest profile impact differences, GSI detects CTP of <500 hPa 

■ However, Aqua MODIS CTP valid at concurrent time as AIRS observation 
indicates CTP is >800 hPa (right circle) and 950 hPa (left circle) 

T (K) ID at <j = 39 (*500 hPa) for 0000 GSI CTP for 0000 UTC analysis on MODIS CTP valid 2240 UTC on 22 
UTC analysis 22 November 2011 22 November 2011 November 2011 
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■ Limited radiance assimilation around 500 hPa in area of largest profile impact 

■ A number of observations retained in the thinning process are not used in the 
analysis due to CTP in GSI being at a higher elevation 

■ Locations of retrieved L2 profiles are larger in number (no data thinning) but also 
provide more data in regions where CTP is lower than 500 hPa 


T (K) ID at o=39 (~500 hPa) for 0000 Assimilated AIRS Radiance Assimilated L2 Profile Locations 
UTC analysis 22 November 2011 Locations at 722cm 1 (~501 hPa) at 500 hPa 



Temperature Innovations 
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■ Unrealistic innovations not the cause of large analysis impact from the profiles in 
this region 

■ Combination of radiances removed due to cloud check and spatial thinning are 
the likely causes for analysis differences 

■ Further investigation into positive or negative analysis and forecast impact 


T (K) ID at c t= 39 (*500 hPa) tor 0000 AIRS Radiance T (K) Innovations AIRS L2 Profile T (K) Innovations 

UTC analysis on 22 November 2011 (OB-BG) at 722cm 1 (*501 hPa) (OB-BG) at 500 hPa 



Summary/ Future Work 
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Summary 

• Parallel experiments using AIRS LIB and L2 retrieved profiles were run for 29 case study 
days for early Winter 2011 

• Forecasts over and downstream regions of low, opaque cloudy regions yield improved T and 
Z anomaly correlations when non-thinned set of profiles is assimilated instead of radiances 

• Initial results indicate that GSI does a good job on the whole of determining cloud-free 
radiances there are some areas coincident with areas of larger profile impact that are 
misrepresented (compared to MODIS) that may result in reduced analysis impact 

Future Work 

• Investigate regions where AIRS radiances have larger impact for possible cloud 
contamination affects 

• Produce quantitative statistics comparing GSI CTPs with MODIS CTPs 

• "Turn knobs" within GSI to determine analysis/ forecast impact from different cloud 
detection, quality, and spatial thinning options 
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